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Abstract: We initiate the systematic computation of the heat-kernel coefficients 
for Laplacian operators obeying anisotropic dispersion relations in curved space- 
time. Our results correctly reproduce the limit where isotropy is restored and special 
anisotropic cases considered previously in the literature. Subsequently, the heat ker¬ 
nel is used to derive the scalar-induced one-loop effective action and beta functions 
of Horava-Lifshitz gravity. We identify the Gaussian fixed point which is supposed 
to provide the UV completion of the theory. In the present setting, this fixed point 
acts as an infrared attractor for the renormalization group flow of Newton’s constant 
and the high-energy phase of the theory is screened by a Landau pole. We comment 
on the consequences of these findings for the renormalizability of the theory. 
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1 Introduction 

Quantizing classical general relativity results in a quantum field theory that is per- 
turbatively non-renormalizable [1, 2]. This feature can readily be traced back to the 
negative mass-dimension of Newton’s constant which entails that an infinite number 
of counterterms is needed to absorb the divergences in the loop-expansion. In terms 
of Wilson’s modern viewpoint on renormalization this is equivalent to the fact that 
the Gaussian fixed point, representing the free theory, is not an ultraviolet (UV) at¬ 
tractor for the renormalization group (RG) flow of Newton’s constant [3]. Contrary 
to QCD, general relativity is not asymptotically free. 

Horava-Lifshitz (HL) gravity [4] (also see [5-8] for reviews) proposes that the 
UV completion of general relativity is provided by a different fixed point, namely an 
anisotropic Gaussian fixed point (aGFP) [9]. The central idea of this proposal is to 
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treat time and space on a different footing. At high energies, local Lorentz invariance 
is broken by imposing an anisotropic scaling between space and time of the form 

t^b~ z t, x'^b-'x*. (1.1) 

This scaling property allows the inclusion of higher-order spatial derivative terms 
in the action while retaining two time derivatives in the kinetic part. This way 
the theory is manifestly unitary. If the anisotropy parameter z equals the number 
of spatial dimensions d, the mass dimension of Newton’s constant is zero and the 
coupling is power counting marginal. This feature makes HL gravity power-counting 
renormalizable and constitutes an important (though not sufficient) ingredient for the 
theory being asymptotically free at the anisotropic Gaussian fixed point (aGFP). The 
quantum properties of HL gravity beyond power-counting arguments have recently 
been investigated analytically [10-14] as well as within Monte Carlo simulations 
based on Causal Dynamical Triangulations [15-18]. 

Since the anisotropic scaling (1.1) treats space and time in a manifestly different 
way, the natural choice for the fundamental degrees of freedom in HL gravity are the 
component fields appearing in the Arnowitt-Deser-Misner (ADM) [19] decomposition 
of the spacetime metric 1 

ds 2 = N 2 dt 2 + <Jij ( dx 1 + N l dt) ( dx J + N J dt,j . (1.2) 

Here N(t,x l ), N l (t,x l ) and x*) denote the lapse function, the shift vector and 
the spatial metric on d-dimensional spatial slices of constant time t, respectively. 
Their scaling dimension is given by 

[IV] = 0, K-] = 0, [N']=z- 1. (1.3) 

The anisotropic scaling (1.1) implies that the symmetry group of general relativ¬ 
ity, comprising spacetime diffeomorphisms, is reduced to so-called foliation preserving 
diffeomorphisms 

1 1 —^ t(t ), x l H > x l {t , x l ) . (1.4) 

The fundamental action for HL gravity then comprises all power-counting relevant 
and marginal interactions that can be built from the field content and are compatible 
with these symmetries. Schematically, this action has the form 

Shl = -^Q J dtd d xN x fd [KijK* - XK 2 + V] . (1.5) 

where K tJ is the extrinsic curvature (see Sec. 2) and K denotes its trace. The poten¬ 
tial V is constructed from the intrinsic curvature tensors and their spatial derivatives, 
together with spatial derivatives of the component fields in (1.2). 

1 We will work in the Euclidean signature. Since the foliation structure naturally defines a global 

time direction, one can safely perform a Wick rotation to change the signature. 
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The precise form of (1.5) depends on the version of HL gravity under consider¬ 
ation. Projectable HL gravity restricts the lapse function to be independent of the 
spatial coordinates x l in order to match the symmetries (1.4). In the non-projectable 
version the lapse is allowed to depend on both the time and spatial coordinates. The 
“healthy extension” of non-projectable HL gravity [20, 21] constructs the potential V 
from the intrinsic curvature as well as spatial derivatives of the lapse function. Since 
the stability of the theory at the quantum level may require the terms added in the 
“healthy extension” we will also comment on this case in the sequel. 2 

Despite the vast literature on HL gravity, little is known about the quantum 
aspects of the theory. The reason for this lack of understanding lies in the unfa¬ 
vorable interplay between the complicated nature of the theory, including a vast set 
of coupling constants, the presence of the anisotropic scaling and the complicated 
way the propagating degrees of freedom are encoded in the formulation. One way 
to bypass these obstructions and address questions related to quantum corrections 
in a curved spacetime in an efficient way, is provided by covariant heat-kernel tech¬ 
niques in combination with the background held method. The key advantage of this 
approach is its universality: once the heat-kernel coefficients are known they can be 
used to study a multitude of physical systems. In particular it is straightforward to 
compute the one-loop effective action of a theory [24, 25]. 

The heat kernel of a second order Laplace-type differential operator A on a D- 
dimensional compact manifold without boundary possesses the Seeley-deWitt (or 
early-time) expansion [25-27] 



( 1 . 6 ) 


where the heat-kernel coefficients ci 2 k are built from curvature tensors and their 
covariant derivatives. A similar expansion also holds for differential operators of 
fourth order [28, 29]. The analogue of (1.6) for differential operators including higher- 
derivatives and a degenerate leading symbol, which is the mathematical equivalent 
of the differential operators naturally appearing in HL gravity, is largely uncharted 
territory, however. The present work therefore serves the following purposes. Firstly, 
the detailed derivation of the heat-kernel coefficients [9] based on off-diagonal heat- 
kernel techniques [30, 31] is given. The result is extended to include terms with 
spatial derivatives of the lapse function which are absent in the projectable version 
of HL gravity. Subsequently, this heat kernel is used to compute the one-loop effective 
action obtained from integrating out scalar fields with anisotropic scaling relations in 
an arbitrary curved background. The one-loop beta functions obtained in this way 
indeed exhibit an aGFP. In this setting the aGFP acts as an infrared (IR) attractor 
while the high-energy phase of the theory is screened by a Landau pole. In this sense 

2 For a recent interesting analysis including mixed derivative terms see [22, 23]. 
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the scalar-driven gravitational renormalization group flow of HL gravity shares many 
features with the one of QED. In particular, the flow is not asymptotically free. 

The rest of the paper is organized as follows. We start by introducing our model 
in Sect. 2. In order to facilitate readability, the heat-kernel coefficients for differential 
operators with anisotropic scaling are summarized in Sect. 3. Owed to the rather 
technical nature, the details of their derivation are relegated to two appendices: 
App. A contains a general exposition of the off-diagonal heat-kernel techniques which 
underlie the explicit computation of the heat-kernel coefficients in App. B. The one- 
loop effective action together with the beta functions for the gravitational couplings 
are obtained in Sect. 4 and we close with a brief discussion of our findings in Sect. 5. 

2 Anisotropic gravity and matter 

2.1 The model 

In this work we focus on matter induced gravitational beta functions. The matter 
sector of our model consists of a (massive) Lifshitz scalar (LS), with action 

S LS = \ J dtd d xNy/a(f)[A t + (A x ) z + m 2 ]0, (2.1) 

where 

A* = -j=d t N~ l y/ad t , A x = ——=d i (r t3 Ny/ad j . (2.2) 

Aya Ayer 

For isotropic scaling, z — 1, the differential operator 

D 2 = A t + (A x ) 2 (2.3) 

in (2.1) coincides with the Laplace-operator Aconstructed from (1.2) with the 
shift vector set to zero 


A 2 0| z=1 = (A t + A x )<j> = A^cj >. (2.4) 

Hence, for z = 1 the matter sector agrees with the one of a diffeomorphism-invariant 
minimally-coupled massive scalar field in a curved background spacetime. Notice 
also that in the projectable case A x coincides with the intrinsic Laplacian A x = 
—a^VjVj, where V,; is the covariant derivative constructed from the spatial metric 
a ij • 

The full bare action of the theory will then be 

•S'[A, N u a, 0] = S H l[N, N u a] + S ls [N , A*, a, 0], (2.5) 

where the Horava-Lifshitz (HL) action is 

Shl = J dtd d xN^ [KijK ij - AA 2 + V] . (2.6) 
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The extrinsic curvature is defined as 


K tj = (2N)- 1 [fan - ViNj - VjNi ], (2.7) 

its trace being K = a l] K l3 . The potential term V is built out of power-counting 
relevant and marginal operators. In the projectable case, one has, in d = 2,3 

V d — g 0 + giR + g 2 R 2 , (2.8) 


and 


V' 3 — go + g\R + g 2 ,iR 2 + g>:>Ri;R'' — gz,iRA x R — g-s. 2 R.ijA x Ry 

+gz,sR 3 + g3 A RRijR ij + g^R^R* • (2.9) 

All intrinsic curvature monomials such as R , R, t] R l] , ... are built from the induced 
metric a i3 on the d-dimensional spatial slices and indices are raised and lowered with 
that same metric. 

In the non-projectable case the potentials include additional terms, e.g., built 
from spatial derivatives of the lapse function 

a i = ^di N - ( 2 . 10 ) 

In [20] for instance the following terms were added to the potential in order to give 
a “healthy extension” of HL gravity 

AV = Mi cii a 1 + u 2 , i R Vj a 1 + u 2 ,2 a* A x a* + u 3t i(A x R) Di a 1 + m 3j2 a* (A z ) 2 a 1 . (2.11) 

This choice limits the potential to terms which contribute to the propagator around 
a flat Minkowski background. In principle, the “healthy extension” of HL gravity 
should include all interaction monomials built from the component fields in (1.2) 
which are compatible with the symmetries (1.4). In total this setup then includes 
0(100) coupling constants, leading to a highly involved theory. 


2.2 Effective actions from the heat kernel 


Starting from an action functional S[cj)\, the one-loop contribution to the effective 
action (EA) takes the form [24, 27] 3 


r x = 


-Tr log 


' 6 2 S ' 
8(t>8(t>_ 


( 2 . 12 ) 


The one-loop effective action is then given by T = S + Ti. Typically, the functional 
trace (2.12) contains divergences which are absorbed into the counterterms contained 

technically one should divide the Hessian by a suitable mass scale fi so that the argument 
of the logarithm is dimensionless. At this stage, this scale will be left implicit. It will appear as the 
renormalization group scale later on. 
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in S. These divergences are crucial for the renormalization of the theory since they 
encode the beta functions governing the energy dependence of the coupling constants. 
A very convenient way to isolate these divergences in curved spacetime is by means 
of the heat kernel. This approach preserves the covariance of the effective action 
with respect to the underlying symmetry group. We will illustrate this technique in 
the remainder of this subsection before applying it to HL gravity in Sect. 4. 

The functional trace of log O , with O an operator (e.g. the Laplacian operator 
on a D-dimensional Riemannian manifold) can be expressed through the eigenvalues 
A n of O as 


^Tr log O = i^logA, 

n 



1 d 

2 ds 


Co(s) 


s=0 


(2.13) 


The last equality defines the zeta-function of O. The zeta-function can be analytically 
continued to the whole complex plane which gives a way to define a finite EA. This 
is the idea of zeta-function regularization [32], There is however another way of 
expressing the functional trace. The sum over eigenvalues appearing in (2.13) can 
formally be written as 

1 poo _ -j poo 

x »‘ = w~\ dt E e ““" = rn / dt • < 2 - 14 ) 

1 Jo „ 1 Jo 



In this way one obtains the formal equivalence 

1 i r°° rlf 

-Tr log O — — / —Tre- t0 (2.15) 

2 2 Jo 1 

which is exact up to a field-independent infinite constant which does not contribute 
to the EA. The relation (2.15) reduces the computation of the one-loop effective 
action to determining the quantity 


K 0 (s) = Tr e~ s ° . 


(2.16) 


Formally, Kq(s) satisfies the heat equation 

d s K 0 (s) + OK 0 (s) = 0, (2.17) 

with boundary condition 0) = 1. For this reason Kq(,s) is called the heat kernel 
of the operator O. 

There exist powerful methods to compute the heat kernel in terms of invariants 
on the manifold [25-27]. The most important one for computing the effective action 
is the Seeley-deWitt expansion which constitutes an asymptotic expansion of (2.16) 
for s —y 0. For the second order Laplace operator this expansion reads 

Tr e~ sA{D) = (4tt) -t s _ ^ f d D Xyfg ^ s " fl 2n,i ^2 ( 2 - 18 ) 

n>0,i 


6 






where the a, 2 n ,i are .D-independent numerical coefficients and the span a basis of 
curvature monomials built from 2 n covariant derivatives. The diffusion time s has 
mass dimension —2, so each term within the sum is dimensionless and the rhs has 
mass dimension zero. 

When we plug the expansion (2.18) into (2.15), we see that there are UV di¬ 
vergences associated with the lower boundary of the integration interval. These 
divergences arise from monomials which come with a negative or zero power of s. 
For monomials multiplied by positive powers of s the integral vanishes at the lower 
boundary indicating that these terms do not give rise to UV divergences. The for¬ 
mer class of interaction monomials precisely correspond to the relevant (s < 0) and 
marginal (s = 0) operators of the theory which require counterterms provided by the 
bare action. 

For HL gravity, the operator we are interested in is the anisotropic Laplacian 
(2.3) which constitutes the Hessian for the action of the Lifshitz scalar (2.1). Based 
on the Seeley-deWitt expansion for this anisotropic operator, we can immediately 
construct the matter-induced effective action at the one-loop level and obtain the 
matter-induced beta functions for the relevant and marginal coupling constants of 
the theory. 

3 Heat-kernel results for anisotropic differential operators 

Based on dimensional analysis, the Seeley-deWitt expansion for the anisotropic dif¬ 
ferential operator (2.3) must take the form 

Tre _s£)2 = (47 t) s f dtd d xN\/a ^ s 1 ^~ L+ ^ + ^ tr. (3.1) 

l,m,n> 0 

Here the number of shift vectors N l , time derivatives, and spatial derivatives con¬ 
tained in the coefficient h[ m n is given by l, m , and n, respectively and tr denotes a 
trace over internal indices. Comparing eqs. (2.18) and (3.1) one first notices a change 
in the overall s-dependence of the rhs as well as the appearance of the anisotropy 
parameter £ in the factors s multiplying interaction monomials including the lapse 
N l and spatial derivatives. This change of scaling is understood as follows. With 
our scaling conventions, time derivatives get mass dimension z, while spatial ones 
keep their standard dimension 1, so that s has dimension — 2z. Then, as before, the 
overall scaling has to compensate the dimension of the integration measure, while 
the powers of s within the sum compensate those of the operators, so that all terms 
appearing on the rhs are dimensionless. 

In App. B we explicitly compute the coefficients b/ m for all interaction monomi¬ 
als containing two time derivatives, two spatial derivatives as well as for all intrinsic 
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Table 1. Heat-kernel coefficients for the anisotropic differential operator D 2 multiplying 
the curvature invariants given in the first column. The second column indicates that 
the corresponding term scales as s~ dim . The heat-kernel coefficients for interaction terms 
typically considered in the framework of HL gravity in d = 2 and d = 3 are given in columns 
3 to 5 and 6 to 10, respectively. The dash implies that the corresponding monomial does 
not appear in the basis of curvature monomials in the given dimension. The coefficients 
for z = 1 agree with the heat-kernel coefficients of the isotropic Laplace operator of second 
order on a manifold with dimension d = 2 and d = 3. 


curvature monomials TZ^n constructed from the metric on the spatial slices a iy Re¬ 
stricted to this (rather general) subclass of terms, the general formula (3.1) becomes 

Tr e~ sD ~ (47t)'~s—“ / dtd d xN^x 

\s n .1 (3.2) 

- (ei K 2 + e 2 KijK ' 3 ) + s z b n a 2n ,i + s z c± did 1 . 

n> 0 i 


Here the a 2ni j are the standard Seeley-deWitt coefficients associated with the intrinsic 



























































curvature monomials TZ^n in the expansion (2.18) (see [33] for an explicit list). The 
coefficients e, and b n depend on d and z and encode the modifications induced by the 
anisotropic scaling of the differential operator. We now discuss the three sectors ap¬ 
pearing in (3.2). The coefficients most relevant for HL gravity are collected in Tab. I. 4 


Extrinsic curvatures 

The coefficients e\ and e 2 associated with the kinetic terms built from the extrinsic 
curvature (2.7) are computed in App. B.2 and read 


e\(d,z) = 


d - z + 3 T(^) 

d + 2 zvU) 


e 2 (d, z) = 


d + 2z T(^) 
' d + 2 zT{%) 


(3.3) 


For z — 1 these equations reduce to ej | 2=] = 1 and e 2 | 2 =i = — 1. This agrees with 
the heat-kernel coefficient a 2 of A ^ expressed in terms of intrinsic and extrinsic 
curvatures and provides a non-trivial check of our result. At criticality, d = z, the 
coefficients (3.3) satisfy the relation 


1 

^1 \z=d "7 ^2 z=d • 

a 


(3.4) 


As we will see in the next section, it is this relation which underlies the fixed point 
structure of HL gravity. 


Intrinsic curvatures 

The coefficients b n are defined through the integral representation (B.46) derived in 
App. B.3. This integral can be evaluated via Mcllin transform techniques giving 

b n(d, z) = l [ d l z 2n +l) , 0 < n < [d/ 2J , (3.5) 

1 ^ 2 +1 ' 

b n(d,z ) = r( ff n fc +fc) / 0 °° dx x d / 2 ~ n+k ~ x (, d x ) k e~ xZ , n > [d/2\ . (3.6) 

Here k = n+1— [d/ 2J and the symbol denotes the floor function, namely 
the largest integer not greater than x. For d even the integral in eq. (3.6) can be 
evaluated explicitly, yielding 


b n (d,z) = (-1 }n-d/2+n/z-d/(2z) 


F (2 n=d + 

■W+7 


n/z—d/(2z),q • 


q E N. 


(3.7) 


Note that the Kronecker-delta puts severe constraints on the intrinsic curvature terms 
that appear in the heat-kernel expansion. 

4 The coefficients a 2n stated in the third and sixth column of the table explicitly depend on the 
dimension d. This is owed to the fact that the results are given for “manifolds of low dimension”, 
where the general basis of geometrical invariants is degenerate. Eliminating the redundant basis 
elements by imposing the vanishing of the Weyl tensor (d = 3) and the vanishing of the Weyl tensor 
together with the Einstein condition (d = 2) generates this d-dependence which is absent in the 
general heat-kernel expansion. 
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At this stage the following comment on the convergence of the integral (3.6) is 
in order. Due to the exponential suppression of the integrand at the upper boundary 
of the integration interval, possible divergences can arise at the lower boundary 
only. For z being integer, the integrand is given by a finite polynomial times the 
exponential. The leading contribution for x —> 0 is proportional to £ d / 2 -L rf / 2 J which 
is manifestly non-negative. Thus for z G {M/0} the integral converges and the b n (d, z ) 
are finite. When z > 0 is non-integer, however, the leading contribution for x —> 0 
can acquire an exponent which is smaller or equal to —1. Acting with k derivatives 
on the exponential, the term which dominates for x —> 0 is of the form 

,1k ( k ~ l \ 

^‘ = - (II - OJ + ■ ■ ■ , (3.8) 

where the dots stand for all terms with a higher power of x. Thus (3.6) becomes 

= fiSi) (iD - o) /“ * +...), (3.9) 

where again the dots denote terms which are subleading at the lower integration 
boundary. The integral is convergent as long as d/2 — [d/2\ — k + 1 + z > 0, or, 
using the definition of k, d/2 — n + z > 0. This establishes that for non-integer z the 
heat-kernel coefficients b n (d, z ) for n < d/2 + z are finite. For n > d/2 + z , the b n (d, z) 
diverge, however. This demonstrates that the operator D 2 with non-integer values 
z does not possess a well-defined Seeley-deWitt expansion due to the divergence of 
the coefficients b n (d, z) with n > d/2 + z. 


Terms containing spatial derivatives of the lapse function 

The coefficient c\ is computed in App. B.4 and reads 


Ci(d, z) = 


i3 z —ir(^ + i) 
e d r (i) 


(3.10) 


Notably, C\ is proportional to z — 1 and thus vanishes in the isotropic case. This is 
consistent with the fact that for z —1 the operator D 2 coincides with the standard 
Laplace operator whose heat-kernel expansion can not contain terms which break 
the invariance under diffeomorphisms. At criticality C\ reduces to 


I 13 (i - 1 r (g - 1) 

1, " d 6 d r(|) 


(3.11) 


The computation of cq serves as an illustration that the off-diagonal heat-kernel 
methods developed in this work are also capable of determining heat-kernel coeffi¬ 
cients relevant for the “healthy extension” of HL gravity. Comparing the complexity 
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of the derivations in App. B.2 and B.4 reveals, however, that finding the heat-kernel 
coefficients for non-projectable HL gravity is computationally much more demanding 
than in the projectable case. Determining all heat-kernel coefficients relevant for the 
“healthy extension” will require a significant computational effort which is beyond 
the scope of the present paper. 

We close this section with the following remark. We explicitly verified that the 
isotropic limit z = 1 of our heat-kernel coefficients correctly reproduces the expansion 
(2.18) written in terms of ADM fields. Moreover, our formulas agree with special 
cases for z — 2 and d = 1 studied in [34], the results for the conformal anomaly 
in z = 2,d — 2 obtained in [35] and the scaling analysis for anisotropic Laplace 
operators in flat space [36]. This match provides a highly non-trivial test for the 
results obtained in this work. 


4 One-loop effective action and beta functions 

We now use the Seeley-deWitt expansion (3.2) to construct the one-loop effective 
action and beta functions for the gravitational sector of HL gravity resulting from 
integrating out the anisotropic scalar field (2.1). 


4.1 General beta functions 

Applying the general formulas (2.12) and (2.15) to the scalar action (2.1), the result¬ 
ing one-loop contribution to the effective action is 


r, 



e- sD2 ' 

Jo s 

- 


(4.1) 


Evaluating the operator trace via the Seeley-deWitt expansion (3.2), one obtains 


r, = - 


(4vr)' 


< 2+1 


dtd d xNy/a 


ds 


2 l+<2/z 

e - sm 2 s -_U_ 


X 


I ( ei K 2 + e 2 ^ s n ' z b n Y, 02 nMl + s* 


Cl a*a + ... 


n> o 


(4.2) 


where the dots indicate further terms containing spatial derivatives of the lapse 
function whose heat-kernel coefficients are unknown. UV divergences are linked to 
the lower boundary of the s-integration. In order to make them visible, we introduce 
a UV cutoff A, so the lower limit is now A~ 2 , and integrate from the cutoff scale A 
down to some renormalization scale /i. 
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It is now straightforward to determine the divergent parts of Id. In the massless 
case, m = 0 these are given by 


r i,A = - 


(47t) 2 

2 

[{d+z)/2\ 

E 


dtd d xNy/a 

2z 




-{ ei K 2 + e 2 K l] K+ 


n= 0 


d — 2n + z 


d—2n+z d—2n+z 

A — fi z 


bn ^ ^ ^ 


W 

2 n 


d — 2+z d — 2-[-2 . 

A * — fi * 1 ci Gqa + ... 


d — 2 + ^ 
subject to the substitutions 

(i - !) 1 ( A " _1 '->■ log(A//x) , if d = z, 


, x 1 / . d— 2n+z d —2n+z 

^ (d - 2n + ^) _1 (A ~ 


H y log (A//i) , if d + 2 = 2n 


(4.3) 


(4.4) 


applicable when a term in the heat-kernel expansion is independent of s or, equiv¬ 
alently, if the interaction monomial is power-counting marginal. In the massive 
case, m ^ 0, the s-integrals give rise to incomplete Gamma functions, T(a, &) = 
f b °° dt f a_1 e“h The analogue of (4.3) then reads 


IYa = 


/ . \_ d+l 

(47t) 2 

2 „ 
L(d+z)/2j 

E Y r < 

n=0 


dtd d xN yfa 


\ (r(i=^i, $) - r(i^4, )) (e, A' 2 + e 2 A'yA'«) 


2 n—d—z m 2 
2z ' A 2 




2 n—d—z m 2 
2z ‘{i? 


bn ^ ^ ® 1 2n,i ^ 


(0 

2n 


+ z r <T. fn) - Ut. 5r) Cl + '' 


(4.5) 


In the next step the UV divergences are reabsorbed by the counterterms present 
in the bare action (2.5). The (scale-dependent) renormalized couplings for the mas¬ 
sive case can then be read off from (4.5). For the couplings appearing in projectable 
HL gravity 


1 

167 tG r (h) 
^r{ aO 

16vrG R (/i) 

9n,i,n(d) 

16nG R (n) 


^2 th / 1—dfz rn 2 ' 


167 tG (47 t ) G + 1 )/ 2 12 

A 


r(- 


167tG 

9n,i 


_i_ ^p/ 1 -d/z m“ \ 

(4 7 r)( rf + 1 )/2 12 1 2 ’i* 2h 


+ 


167 tG (47t)G+ 1 )/ 2 


% bn @>2n,i k( 


2 n—d—z m 2 ' 
2 z LI 2 2 


(4.6) 


Here G, A, and g V)l are the values of the renormalized couplings at fj, = 0, which 
constitute free parameters of the theory. The renormalized coupling constants in the 
“healthy extension’’ follow the same pattern as the g n ,i,R(l^)- I 11 particular, 


U lAd) U 1 , 1 „ v( 2-d-z rn? 

167T Gr{h) 167 tG (47t)G +1 )/ 2 1 2 Z 1 fi 2 


(4.7) 
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The scale-dependent renormalized coupling constants of the massless case are found 
in an analogous way either by reading them off directly from (4.3) or taking the 
limit m —> 0 in eqs. (4.6) and (4.7). In this limit the incomplete Gamma-functions 
simplify to a power-law running in /i if the corresponding coupling is relevant and 
a logarithmic /i-dependence for marginal couplings. Since the subsequent discussion 
involves only renormalized couplings, we will from now on drop the subscript R for 
ease of notation. 

The scale-dependence of the renormalized coupling constants is conveniently en¬ 
coded in the beta functions of a theory. For this purpose, we use the renormalization 
group scale fi to construct the dimensionless couplings 

AKaO = ^(aO A* ^ j 9n,i{f^) = 9n,i (AO A 4 ^ j U n ,i {AO = ^n,i(/^) A 4 ^ j (4-8) 


where the mass-dimensions are determined through 




(4.9) 


The beta functions /id^g = /3g are obtained by re-expressing eqs. (4.6) and (4.7) in 
terms of the dimensionless couplings and taking the logarithmic /i-derivative. Ex¬ 
plicitly, 


& = 2 >) 9 - | 9 2 , 

h = -| (4x)- ( ' , - 1,/2 e-'‘ _w (e, + Ae 2 ) g. 

Pto'i = - ^in9n,i + ^(47r) _(<i_1)/2 e _M 2m “ (4 b n a 2n ,i — § 9„,ie 2 ) , 
Ah = - 2 & fti + g (4vr) _(d_1)/2 e _M 2m2 (4 Ci - § iq e 2 ) • 


The beta functions for the massless case follow straightforwardly by setting m = 0. 
In this case the beta functions form an autonomous system of differential equations, 
independent of fi. 

The system of equations (4.10) is the central result of this subsection. Owed to 
the fact that it captures the quantum effects arising from integrating out a minimally 
coupled scalar field, the system of flow equations (4.10) is structurally relatively 
simple. The beta function for the dimensionless Newton’s constant is independent 
of all other couplings and can thus be solved on its own. Subsequently, this solution 
can be plugged into the remaining equations which can again be solved individually. 
Thus considering the “healthy extension” does not modify the renormalization group 
flow in the “projectable” sector of the theory. In this sense the flows found for the 
“healthy extension” form a direct extension of the ones obtained in the projectable 
theory. 

We close this subsection with the following remark. Albeit the full set of heat- 
kernel coefficients for the relevant and marginal interaction terms of the “healthy 
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extension” is not known, it is clear that the beta functions for the couplings u n ^ 
will have the same structure as the ones for the couplings g n ^ with different coef¬ 
ficients b n a,2n,i- Parameterizing the unknown coefficients by c n ,* where 2 n is again 
the number of spatial derivatives in the interaction monomials and i enumerates the 
corresponding basis monomials, one concludes that 

Pu nti = -2 u n ,i + g (4vr) _( ^ _1)/2 e"“ m2 (4 c nti - § u nti e 2 ) • (4.11) 

As we will see in the next subsection, this parametric form of the beta functions 
already allows to draw conclusions about the structure of the renormalization group 
flow. 

4.2 Renormalization group flows at criticality 

In terms of HL gravity, the most interesting situation is when the system is at crit¬ 
icality with d = z = 3. In this case, Newton’s constant is marginal and the theory 
is power-counting renormalizable. We will discuss the renormalization group flows 
entailed by (4.10) and the extension (4.11) in this subsection. 

Massless flows 

Restoring the number of scalar fields n s , the scalar-induced massless beta functions 
for the kinetic terms at criticality simplify to 

= A = ^(3A-l )g. (4.12) 

This system is supplemented by the beta functions for the potential couplings in the 
projectable sector 

/%, = _ 2g 0 + n s f (^o + | go) 

P~91 = “ 101 + n sn ( a 2,i + 19i) 13 s 

P§2,i = - f 02,i + n 8 f ( b 2 a 4 ,i + l <72,i) , 

P~93,i =J *4( i 3 «6 ,i + l 93,i) , 
and the “healthy extension” 

Au = - | ui + n s f ( ci + | ui) , 

Pu 2 ,i = - f u 2 ,i + n s f ( C 2 ,i + | W 2 ,i) , ( 4 . 14 ) 

^«3,i = n s^{ C 3 + i « 3ii ) . 

The beta functions (4.12) and the marginal couplings in the projectable sector agree 
with the ones obtained from the functional renormalization group [9]. This match 
is in agreement with the expectation that the one-loop beta functions for marginal 
couplings are universal in the sense that they do not depend on the regularization 
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Figure 1. RG flow in the g—X plane for the massless scalar case in d = z = 3. The vertical 
line represents the isotropic plane A = 1. The arrows point towards the IR, and the black 
dot indicates the position of the aGFP. We clearly see that the aGFP is an IR attractor 
for g > 0 and a UV attractor for g < 0. 

scheme. Typically, this universality does not extent to relevant coupling constants. 
On this basis, it is not expected that the beta functions for g 0 , gi and g^.i agree in 
the two computations. 

The most prominent feature of the system (4.12), (4.13) and (4.14) is the pres¬ 
ence of an anisotropic Gaussian fixed point (aGFP) where all beta functions vanish 
simultaneously 



(4.15) 


a; - «), - o. fi3,i = -5c 3 , ( . 


This fixed point is two-fold degenerate in the sense that it is located at the intersection 
of a line of Gaussian fixed points with g = 0 and a line of non-Gaussian fixed 
points parameterized by the anisotropy parameter z [9]. The aGFP corresponds to 
a free theory where the propagators exhibit anisotropic scaling with z = 3. Notably, 
the position of the fixed point is independent of n s . This implies that the fixed 
point solution becomes exact in the large-n s limit, where the contribution of the 
gravitational loops to the beta functions is negligible. It is actually this fixed point 
which underlies the expected renormalizability of HL gravity both in the projectable 
as well as in the non-projectable version. 

At this stage, it is instructive to investigate the subsystem (4.12) in order to 
understand the properties of the aGFP in more detail. Fig. 1 shows a set of sample 
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1.0 


0.5 



Figure 2. RG flow for g and A as a function of t = In (fi/m). The dashed line represents 
the isotropic plane A = 1. The flow of g and A freezes out for g 2 < m 2 . 


solutions illustrating the renormalization group flow in the g —A plane with the aGFP 
marked by a black dot. This establishes that the aGFP actually acts as an IR 
attractor for trajectories starting with a positive Newton’s constant. The analytic 
solution of the system (4.12) is 


M = 


5x 


5nai + log/i ’ 


A(/i) = 


O 2 


+ 


log fJ, 


5x0,1 + log R 3(57rai + log g) ’ 


(4.16) 


where ai and 02 are integration constants. Thus g increases with increasing en¬ 
ergy scale, and the system becomes strongly coupled at high energy. In fact, the 
dimensionless Newton’s constant apparently becomes infinite at some finite energy, 
resulting in a Landau pole. In this respect, the renormalization group flow of the 
system resembles the features known from QED. 


Massive flows 

The beta functions incorporating a non-zero mass can be obtained from the system 
(4.12), (4.13), and (4.14) through the formal substitution n s 1 —» n s e“ m2 ^ 2 . The 
presence of the mass scale motivates introducing the “renormalization group time” 
t = In {fi/m). Fig. 2 then illustrates the typical renormalization group flow of g(t) 
and A(£). In this case the flow of the coupling constants freezes for g < m. This is 
the standard decoupling phenomenon of a gapped theory. Due to the decoupling of 
the scalars below the scale m 2 , the IR limit of the matter-induced flow is trivial and 
the infrared value of the coupling constants can be tuned to any value by choosing 
suitable initial conditions. An interesting observation is that for energies g > m, the 
flow is still attracted towards the aGFP before the couplings freeze out. This means 
that if the isotropic system is perturbed by a Lorentz breaking (scalar) operator at 
some high energy, the renormalization group flow will drive the system away from 
the isotropic plane. It would be very interesting to understand the phenomenological 
consequences of this instability in more detail. 
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5 Conclusions 


One of the primary motivations for ffL gravity is the prospect of obtaining a quantum 
theory for gravity which is both perturbatively renormalizable and unitary. Its core 
feature is the anisotropic scaling of spatial and time directions (1.1) which, at critical¬ 
ity, leads to a dimensionless Newton’s constant without introducing more than two 
time derivatives in the fundamental action. This renders the theory power-counting 
renormalizable without introducing negative norm states plaguing higher-derivative 
gravity. An open challenge in the program consists of going beyond the power¬ 
counting arguments and establishing that HL gravity indeed possesses a suitable 
renormalization group fixed point which renders the theory asymptotically free. 

In this paper we have devised a general method to derive renormalization group 
flows in HL gravity based on the effective action. Following a strategy developed for 
quantum field theory in curved spacetime (see, e.g., [24, 27]), we used the Seeley- 
deWitt expansion of the heat kernel for anisotropic differential operators to compute 
the one-loop effective action originating from integrating out a (massive) scalar field 
with anisotropic dispersion relation minimally coupled to HL gravity. Our approach 
naturally links with the background field method in curved space and directly leads 
to a covariant effective action. Our technique is very general and can systematically 
be generalized to other types of matter fields including fermions and gauge fields as 
well as to graviton loops. 

The main technical achievement of our work is the first systematic computation 
of the Seeley-deWitt expansion for an anisotropic Laplace-type operator D 2 = A t + 
(A x ) z in a curved spacetime. The heat-kernel coefficients for the two-derivative terms 
and all intrinsic curvature monomials are obtained explicitly for arbitrary dimension 
of the spatial slices d and anisotropy z. Our derivation is based on off-diagonal heat- 
kernel techniques. Notably, the applicability of this method is not limited to a specific 
version of HL gravity: we explicitly demonstrated that the heat-kernel coefficients 
associated with operators appearing in the potential of the “healthy extension” can 
be obtained as well, even though their computation is much more demanding than 
finding the corresponding coefficients appearing in the projectable version of the 
theory. In fact, obtaining the complete set of heat-kernel coefficients associated 
with the power-counting relevant and marginal interaction terms in the “healthy 
extension” may be as complex as a two-loop computation in general relativity and 
will thus requires a dedicated effort. We hope to come back to this point in the 
future. 

The key features of the scalar-induced renormalization group flow at criticality 
z = d = 3 are summarized as follows. In the case where the scalar field is massless, 
the flow exhibits an anisotropic Gaussian fixed point (aGFP) [9]. For the relevant 
couplings of the theory, this fixed point acts as a UV attractor, while it is UV repulsive 
for a (dimensionless) positive Newton’s constant. Tracing the renormalization group 
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flow emanating from the aGFP in the IR towards the UV, the renormalization group 
trajectories terminate in a Landau pole. In the case where the scalar field has non¬ 
zero mass m the renormalization group flow freezes at scales /i 2 < m 2 and the value 
of the coupling constants in the IR constitute free parameters of the theory. Notably, 
these features are independent of the version of HL gravity imder consideration and 
apply to the projectablc and non-projectable versions of the theory. This universality 
is easily traced back to the observation that the different versions of HL gravity 
differ in their gravitational sectors only. Since the present model does not include 
the contribution of graviton loops to the beta functions, our results are independent 
of the precise structure of these sectors. 

The gravitational renormalization group flow induced by the massless scalar field 
is strikingly similar to the one encountered in QED. In both cases the fixed point 
corresponding to the free theory is an IR attractor and the high-energy behavior 
of the theory is screened by a Landau pole. Notably, this behavior becomes exact 
in the large iV-limit when the number of scalar matter fields dominates over the 
gravitational contribution to the beta functions. This does not imply, however, that 
HL gravity cannot be asymptotically free. Similarly to the one-loop beta function of 
QCD with rif flavors, 



(5.1) 


the sign of the beta function in HL gravity may depend on the type and precise 
number of matter fields appearing in the theory. It would be very interesting to 
determine this dependence in detail. 
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A Off-diagonal heat-kernel techniques 


The computation of the Seeley-DeWitt expansion for the operator D 2 crucially hinges 
on the ability to evaluate operator traces of the form 


Tr, 




O 


= (xA e 




o \ Xi ) 


(A.l) 


where = —g^ u D^D u is the standard second order Laplacian on a D — d + 1- 
dimensional manifold, O is an operator insertion typically including covariant deriva¬ 
tives Dfj, constructed from g IJV , and the \xA) form a basis of position eigenstates car¬ 
rying internal indices i. In the present computation, the Laplacian acts on scalar 
fields so that the internal indices are absent. 

The off-diagonal heat-kernel allows to evaluate traces of the form (A.l) as follows 
[30, 31]. 5 Inserting a complete set of states, eq. (A.l) can be written as 

{x\ Oe~ sA(D) |x) = (x| O la/^x'le -8 ^ 1 x) 

f „ (A. 2) 

= J d D Xy/gtr[OH(x,x'-,s)] x=x , , 


where H(x,x'\s) = (x'|e sAtD> |x) is the heat kernel at non-coincident points. This 
“off-diagonal” heat-kernel obeys the asymptotic expansion 

, oo 

H(x,x';s) ~ ( 47 rs) _D//2 e 2a ) ^ s n A 2n (x, x'). (A.3) 

n =0 

Here, a(x,x') denotes half the squared geodesic distance between x and x' satisfying 

= a(x,x) = a ;i j.(x,x) = 0, cr.^x, x) = g^x) , (A.4) 

and A2n(x,x'), n G N, are the off-diagonal heat-kernel coefficients. The latter are 
purely geometrical two-point objects, formally independent of the dimension D. At 
coinciding points A 2n (x,x) agrees with the standard heat-kernel coefficients (B.44). 
Their values at non-coincident points can be encoded in a Taylor expansion around 
x with coefficients given in [38] and references therein. 

In practice, eq. (A.2) is evaluated by acting with the covariant derivatives con¬ 
tained in O on H(x,x']s) and subsequently taking the coincidence limit. In order 
to simplify the occurring structures one can assume that the covariant derivatives 
within O are symmetrized, i.e., 

n 

0 = (A.5) 

k =0 

5 For some early work on the spin-dependence of the off-diagonal heat-kernel also see [37]. 
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with totally symmetric matrices M ai '" a2k . This form can always be achieved by 
successively writing non-symmetric combinations of the derivatives as a sum of sym¬ 
metric and antisymmetric pieces and expressing the latter in terms of curvature 
tensors using commutator relations of the form 


[ -D/x , D v ] CL\ ®<r • 


(A.6) 


Based on the structure (A.5), it then suffices to have an explicit form for the 
quantities 

Hoti...a2n = ^(ai ’ ’ ’ ^a 2 n) H(x, J (A.7) 

where the overline denotes the coincidence limit x —> x' and the round brackets 
indicate the synunetrization of indices with unit strength, (a/3) = \(a/3 + /3a). In 
terms of the 77-tensors the trace (A.l) takes the rather simple form 


Tr 


-sA( D ) 


o 


dPxyfg M ai ' a2k H ( 


a\...a 2 k 


(A.8) 


fc >0 

The explicit expressions for the 77-tensors can be constructed from the covariant 
derivatives of a(x,x') and A 2n (x, x') evaluated in the coincidence limit. Based on 
the properties (A.4) one can prove recursively that a(x, x') satisfies [30] 


cr(x,x').^ ai ... an) = 0, n> 2. 


(A.9) 


Furthermore, a recursion relation based on the heat equation gives [30] 


A 0 (x,x') = 1, A 0 (x,x')^ = 0, A 0 (x,x , ).( Jlu ) = , 

A 2 (x,x') = ± R. 


(A.10) 


Coefficients A 2n (x,x') with n> 2 and terms based on A 0 (x,x') and A 2 (x,x') involv¬ 
ing additional covariant derivatives contain more than two covariant derivatives and 
do not contribute to the computations carried out in this work. They can be found 
in [33]. Taking covariant derivatives of (A. 3), evoking the coincidence limit, and 
substituting the relations (A. 4), (A.9) and (A. 10) yields the following asymptotic 
expansion of the 77-tensors 


H - (SysTi l 1 + ¥ R ) • 

Ha — 0 , 

— (4 7rs )£>/2 { - 2s ^ + + 6^ a/3 } ’ 

H-apiJ, — 0 ? 


Ha ^ v ~ (4 tt s ) d / 2 l 4 s 2 


\ , 2 [QolhQPv 9au9(3fi “I - 9oip9nv) (l “1“ gS-R) 


^2 S ^9oliiR(3v 9cx.v'R(3^l "I - 9f$ii^ j OLV 9pvRa.ii “I - 9apRfiv 9fivRaft) [ 5 


(A. 11 ) 
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which is exact up to terms including three or more derivatives. Eq. (A.3) in combi¬ 
nation with the //-tensors (A. 11) constitutes the central result of this appendix. In 
App. B we use these formulas to evaluate the operator traces encoding the Seeley- 
deWitt coefficients of D 2 . 


B Anisotropic heat-kernel coefficients: derivation 

This appendix provides the technical details underlying the derivation of the heat- 
kernel coefficients for the anisotropic differential operator (2.3) reported in Sect. 3. 


B.l General setup 

The goal is to determine the coefficients appearing in the Seeley-deWitt expansion 
(3.1) of 

K D 2 (s) = Tre~ sD2 (B.l) 

for the anisotropic differential operator D 2 = A t + (A x ) 2 defined in eq. (2.3). A key 
advantage of the Seeley-deWitt expansion of D 2 is its covariance with respect to the 
symmetry group (1.4). The structure of the interaction monomials is then largely 
determined by these symmetries. This feature allows to simplify the computation of 
the heat-kernel coefficients by working with a specific subclass of metrics in order to 
simplify the computational effort. This class has to be sufficiently general in order to 
uniquely identify the curvature monomials of interest. Covariance of the heat-kernel 
then implies that the coefficients determined from this specific set of backgrounds 
carry over to the general case. 

In practice the heat-kernel coefficients of D 2 will be obtained by recasting the 
operator trace (B.l) into the form 


Tr e 


-sD 2 


Tr 


e -sA( D ) 



(B-2) 


where A^ = D 2 \ z= ± is the standard Laplace operator (2.4) and the operator in¬ 
sertions O are built from commutators of A t and A x . The rhs of this expression 
can be evaluated by employing the off-diagonal heat-kernel techniques described in 
App. A. The relation (B.2) is systematic in the sense that at a given order of the 
Seeley-deWitt expansion multi-commutators contribute up to a finite order only. 6 

The first step towards obtaining (B.2) uses the Zassenhaus or inverse Campbell- 
Baker-Hausdorff formula 

e -s(A+B) = e -sA e ~sB e -^lA,B] e -^([A,[A,B]]+2[B,lA,B}]) . . . (g.3) 

6 In the sequel, we restrict ourselves to multi-commutators with at most two spatial or time 
derivatives acting on the component fields. Commutators containing more derivatives of the com¬ 
ponent fields do not contribute to the heat-kernel coefficients computed in this paper and will not 
be considered. 
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to split the exponential 


Tr e~ sE>2 = Tr [e~ sAt e~ s{Ax)z C(A t , A x )] . 


(B.4) 


The operator insertion C has an analytic expansion in s = 0 whose coefficients can be 
found by expanding the commutator terms in eq. (B.3). Retaining all commutators 
which contain up to two time derivatives acting on the component fields 


C(A t , A x ) ~ 1 - j [A*, (A x ) 2 ] - ^ ([A t , [A^ (A*)*]] + 2[(A X )*, [At, (A x )*]]) 
+ ^[A t ,(A x )*] 2 + 0(a 5 ). 


(B.5) 


Generically, these commutators involve non-integer powers of A x . By induction one 
can prove that these are related to the standard commutators of A t and A x by 

[ At, (A x )‘ ] = * (A x ) z ~ l [ A*, A, x ] + \z(z - 1) ( A x ) z ~ 2 [[ A t , A x ], A x ] +... , (B.6) 

where the dots contain multi-commutators of higher order. Without loss of generality, 
we assume that all terms in C(A t ,A x ) are “normal ordered” in the sense that all 
operators A x are at the left and the A t appear at the right, 


C(A,, A.) = v f„(A x )C,ft,(A t ). 


(B.7) 


The Ci contain all parts of the commutators which do not combine into the Lapla- 
cians. Using the cyclicity of the trace, eq. (B.4) is cast into the form 


Tr e ~ sD2 = Tr "" S(A * )2 


f Xti (A x )C 


(B.8) 


The goal of the next step is to eliminate the exponent z appearing in 
For this purpose, we introduce the anti-Laplace transform W(v) of a function W(x 
as 

/»oo 

W{x)= / dvW{y)e~ xv . 


Setting 


one arrives at 


WAx) = e 


Tr e" 


E 


fx,i ( X ) 

) 

dvW,(v) Tr / W (A.) e~‘ A ‘ e~" A ’ C 


(B.9) 

(B.10) 

(B.11) 


In principle, the functions may also be eliminated by a suitable Laplace transform. 
Since f t . t = ( A t ) ni with n* being integer, it is more convenient to generate these terms 
by taking suitable derivatives of the trace with respect to s 


Tr e 


-sD 2 


'LL dvW ‘( v A~T s ) ^ 


d 


e~ sAt e~ vAx C 


(B.12) 
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In the next step, one introduces an auxiliary Laplace operator 

A X = ~A X . (B.13) 

s 

This auxiliary operator is associated with the auxiliary spatial metric 

s 

a ij = -a ij . (B.14) 

Writing (B.12) in terms of this auxiliary Laplacian allows to combine the two expo¬ 
nentials by applying the Campbell-Baker-Hausdorff formula in its standard form 

e -sA te -sA x = e-dAt+A*)^^ (B.15) 


with 


2 3 

I3(A t , A x ) ~ 1 + — [A t , A x ] + — ([A t , [At, A x ]] + 2[A X , [At, A a 


■[At, A x ] 2 + (9(s 5 


(B.16) 


Substituting this result into (B.12), the rhs indeed assumes the form (B.2) with the 
D-dimensional Laplacian = A t +A x defined with respect to the auxiliary spatial 
metric (B.14) 


Tr e 


-sD 2 


£ 


dvWAv) 




(B.17) 


This operator trace is then readily evaluated by applying the off-diagonal heat-kernel 
techniques of App. A. Expanding all curvatures into d-dimensional ADM invariants 
by using the Gauss-Codazzi relations and converting back to the original metric a l3 
the rhs is readily expressed in terms of curvature monomials. 

The final u-integral can be evaluated using Mcllin-transform techniques, see, e.g. 
[39, 40]. Using 

/•OO 

Q n [W}= / dvv~ n W(v), 

Jo 

dxx n ~ l W(x ), 

the integrals over the anti-Laplace transform W{y) are related to moments of the 
initial function W(x). Combining eqs. (B.5) and (B.6), one finds that the func¬ 
tions Wi{x), eq. (B.10), fall into the classes collected in Tab. 2. Their Q-functionals 
coincide with integral representations of the Gamma-function, which allows deter¬ 
mining the heat-kernel coefficients in closed form. 



n > 0 

(B.18) 
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Wi(z) = e~ sxZ , 

W 2 (x ) = szx^e-™*, 

W 3 (x) = s 2 z 2 a; 2 ^ 2 e _S3;Z , 
W 4 (x) = ^ sz(z — 1) x z ~ 2 e~ sx 


fo°° duu d/2 Wd(u) = s d/(2z) 

J 0 °° dv V ^ d/2+V> W 2 (v) = S~ d/( - 2z) T ^!{d/2) 

f°° d vv -W 2 +2) - s -d/(2*) r W( 2 d) rf±2£ 

JO UL L VV3\U) * zT(d/2) d+2 

r°°rfwu- (d/2+2) W^lV) = s - rf /( 22 ) r (° ! /( 2 ^) £ni 

Jo UU U vv 4\ u ) b zT{d/2) d+ 2 


Table 2. Integrals occurring in the evaluation of the traces (B.24) for the extrinsic cur¬ 
vatures and monomials of the healthy extension. The functions Wj are related to the 
functions Wj via the integral transform (B.9). 


B.2 Kinetic (extrinsic) terms 

In order to obtain the heat-kernel coefficients built from the extrinsic curvature (2.7) 
it is convenient to choose a class of backgrounds A4 = S 1 x t T d where the radii of 
the tori depend on t. Explicitly, we set 

N = 1, N l = 0 , (Tj j(t, x) = (Tij(t ). (B.19) 


In this case the intrinsic curvatures and the derivatives of the lapse function vanish 
while the extrinsic curvature simplifies to K i3 = The operators (2.2) take the 

form 

A t = - a^Kijdt , A, = -a ij 8,8 3 , (B.20) 

while the non-vanishing Christoffel symbols of the the D-metric on A4 are given by 

r%- = —Kij, r tj = <j ik K k j . (B.21) 


Evaluating the trace (B.l) on the background (B.19) projects the rhs of the master 
equation (3.1) onto the terms constructed from K 7J 


Tr e 


-sD 2 


M 


~ (47?) 


~(d+ 1)/2 


dtd d x\/a [b 0 + |ei K t] K lJ + |e 2 A' 2 + ...] 


(B.22) 

where the dots indicate terms containing more than two time derivatives. The aim 
is to determine the unknown coefficients b 0 , e\ and e 2 . 

This computation proceeds along the lines outlined in the previous subsection. 
Starting from (B.20) the commutators entering into (B.6) and (B.16) are given by 


[At, A x ]<f>=- 2(2 K*d t + K ij + KK lJ )d l d j (p , 

[A t , [A t , AJ] </> = 8 k^d^didj ^, (B.23) 

[A x , [A*, A J] </> = 8 K i *K lk d i d j d l d k (l >, 


where the dot indicates a partial derivative with respect to t and all terms with more 
than two derivatives acting on component fields have been dropped. 
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The next step constructs the explicit form of (B.4) by substituting (B.5). Con¬ 
veniently, the result is decomposed according to 


Tr e~ sD2 ~ Ti + T 2 + T 3 + T 4 + T 5 , 


(B.24) 


with one trace associated with each commutator 
Ti =Tr [e~ sAt e~ s ^ Y ] , 

T 2 = _ ^Tr[e- sA ‘e- s ( A ^ [A t , (A *) 2 ]] , 

Ts [ e - A . e -(A.)‘ ([A t , (A-*,.) 2 ]) 2 ] , (B .25) 

Ta=~j Tr[e~ sAt e- s ^ z [(A*) 2 , [A t , (A*) 2 ]] ] , 
Tr[e- sAt e- s ^ z [A t , [A t , (A*) 2 ]] ] . 

The commutators appearing in these expressions are subsequently expressed through 
the basic commutators (B.23) via the relation (B.6). Investigating the resulting 
operator insertions, one finds that the traces T 3 , T 4 , and T 5 already contain the 
maximum number of two time derivatives acting on spatial metrics cr ?J . At this point 
all operators can be commuted freely and the partial derivatives may be promoted 
to covariant derivatives with respect to the D-metric on AT The commutators 
and connection pieces include further powers of the extrinsic curvature and therefore 
do not contribute to the present computation. 

We illustrate these steps by the detailed computation of T 3 . The substitution of 
the explicit commutators gives 

T 3 ~ 2 sV Tr [ e~ sAt e~ s{A * Y (A *) 2 * -2 K ij K kl didfikdt d 2 \ . (B.26) 


Subsequently, the trace is written in terms of the Laplace anti-transform W 3 of 
the function W 3 (x) = s 2 z 2 x 2z ~ 2 e~ sxZ and the partial derivatives are completed into 
covariant derivatives D M 

/»oo 

T 3 ~ -2s 2 / dv W 3 [y) Tr [A t e~ sAt e~ vA * K ij K kl D.D^DM] . (B.27) 

Jo 

The operator insertion A t is then converted into a derivative of the trace with respect 
to s and the two exponentials combine into the auxiliary Laplacian A^ D ^ constructed 
from the rescaled spatial metric (B.14) 


T 3 ~ —2s 2 


dvW 3 (v) 



Tr 


e ~ sA(D) K ij K kl DiDjD k Di 


(B.28) 
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At this stage the trace can easily be evaluated using the off-diagonal heat-kernel. 
Substituting (A. 11 ) with the curvature terms set to zero gives 


T 3 ~ 2s 2 / dvW 3 (v)d s 


'0 


( 47 rs) 


d +1 


dtd d xVaK tJ K kl ^ ( a ik aji + aud jk + d i:j d k i) 

(B.29) 


Converting back to the initial spatial metric cr via (B.14) and taking the s-derivative, 
one arrives at 

/»OO p 


T 3 ~ -i (4 tt)' 


J dv v~^ 2+2) W 3 {v) J dtd d x^ [2K ij Kij + K 2 ] . 

(B.30) 

Finally, the integral is carried out with the help of Tab. 2, yielding the final form of 
the trace 

d 


T 3 ~ 


4 


(4vr)- (d+1)/2 s -l( 1+ci / 2 ) f dtd d xyfcy s [2KijK ij + K 2 ] . 

* T( 2 ) d + 2 J 


(B-31) 


Computing T 4 and I 5 proceeds along the same lines which results in 

d 


T,~-\ (4vr )- (d+1)/2 s ~^ 1+d/z) [ dtd d xy/as [2K^ + K 2 ] , 

3 zl(f) « + 2 J 

( 4 vr )" (d+1)/2 QAl d + 


^r(f) d + 2 


dtd d x^/a s 2I\ijI\' J + K 2 ] . 


(B.32) 


The evaluation of T\ and T 2 is slightly more involved since there are additional 
contributions from the £btype commutators (B.16) and the completion of the partial 
into covariant derivatives. Upon applying (B.15) the trace T\ becomes 

,,2 


Ti ~ / dv Wi Tr e 


-sAW 


1 + ~2 [^t, 


+ - g _ ([At, [At, A x ]] + 2[A X , [A t , A x ]]) + — ([A t , A a 


(B.33) 


with the function Wfi defined in Tab. 2. Substituting the explicit commutators (B.23) 
and applying the off-diagonal heat-kernel, the final result for T\ is 

d \ r 

dtd d xyfo [1 + f (K 2 - K t] K ir )] 

(B.34) 


T ~ f47T1 — (rf+l)/ 2 s ~\0-+d/z) 2z) 

1_ zT{i) 


+ A( 4vr )-U+ 1 )/ 2 s -|(i+^) [ dtd d xy/^s \2K i:j K ij + K 2 }, 

£ T(|) J 

where the first and second line captures the contribution from the “one” and the 
commutator terms, respectively. Finally, the evaluation of T 2 yields 

T 2 ~ \ (4vr )" (d+1)/2 s ~^ 1+d/z) [ dtd d Xyfc? s [2K i:j K ij + K 2 ] . (B.35) 

2 2 1(f) d + 2 J 
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Summing up the contributions (B.31), (B.32), (B.34) and (B.35) gives the final 
result for the heat kernel evaluated on Af 


Tr e 


-sD 2 


M 


K ' znh 


dtd d xy/a [1 + | (K 2 - K tJ K lJ )] 


(4 J r)-' d + 1) '' 2 d+«/=) 11^1 j=i [ dtd d x\/a [| (K 2 + 2K ij K ij )] . 


*r<5) 


(B.36) 


Comparing this result to the initial ansatz (B.22) allows to read of the coefficients 


_ £(j|) ; _ d — z + 3 T(A) _ d + 2z r(^) 

~ zT{%)' &1 ~ d + 2 Z T{±)' 62 “ d+2 z Y{f)' 


(B.37) 


This result completes the derivation of the heat-kernel coefficients appearing in the 
kinetic sector of HL gravity. 


B.3 Potential (intrinsic) terms 

In this subsection we determine the heat-kernel coefficients associated with monomi¬ 
als constructed from the intrinsic curvature. In this case it is convenient to evaluate 
the trace (B.l) on the product space A 4 = S 1 x A4 where At is a d-dimensional 
compact manifold coinciding with the spatial slices. The component fields of this 
background are taken as 


N = 1 , N l = 0 , crij(t,x ) = cTij(x) . (B.38) 


so that all fields are independent of time. The operators (2.2) in this background are 
given by 

A t = -d 2 , A . x = -a ij V i V j . (B.39) 

where V* is the covariant derivative constructed from a l3 , i.e., V^cr^ = 0. 

A direct consequence of the background (B.38) is that the operators (B.39) 
commute 

[ A t , A x ] = 0. (B.40) 

The evaluation of the trace makes essential use of the product structure 


Tie~ sD — Tr [e- sAt l • Tr [e“ s(A:c) H 

M 

POO 

= Ti[e~ sAt ] • / dvW!(v)Tr[e- vA *] . 

Jo 

The first trace is the heat kernel on S 1 , 


(B.41) 


Tr [e- sAt ] = 


(47rs) 1 / 2 


dt, 


(B.42) 
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while the second trace has the standard Seeley-deWitt expansion of the Laplacian 
on the d -dimensional manifold M. 


Tr [ e = * / ddx Y, yna2n • ( R43 ) 

^ J n>0 

The coefficients a, 2 n can readily be found in the literature [25-27, 33]. In particular 
the first two are given by 


<2.0 — 1 j a 2 — • • • j (B.44) 

with R the (intrinsic) curvature on A4. For manifolds of dimension d = 2 and d = 3 
where the curvature tensors obey special relations, the relevant coefficients d 2 n are 
listed in columns 3 and 6 of Tab. 1. 

Substituting (B.42) and (B.43) into (B.41) we find an expansion of the form 


(B.45) 


where the coefficients b n are given by 

/»oo 

b n (d,z) = s ^ / dvv-V^W^v). (B.46) 

Jo 

Here the s-dependent prefactor has been chosen in such a way to cancel the s- 
dependence of the integral, so that the coefficients b n (d, z ) are independent of s. The 
explicit value of the integral may again be found by applying the Mellin transform 
techniques (B.18) which leads to the results (3.5) and (3.6). 


B.4 Non-projectable potential terms 

Finally, we determine the heat-kernel coefficient of the interaction monomials in¬ 
cluding spatial derivatives of the lapse function. Similarly to the case of the intrinsic 
curvatures, we limit ourselves to terms containing up to two derivatives of the com¬ 
ponent fields. These terms can be determined on a background A4 = S 1 x T d with 
the component fields restricted to the form 

N = N(x), IV* =0, ( r ij (jt,x) = 8 ij , (B.47) 

where the lapse N(x) depends on the spatial coordinates x only and the spatial 
metric is given by the Kronecker symbol. Setting a* = d{ In N, the two operators 
(2.2) reduce to 

A, = - N~ 2 d 2 t , A x = -ftdidj -a'di, (B.48) 

and the non-zero Christoffel symbols of the D-dimensional geometry are 

r \t = -N 2 a l , r t it = ca. 
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(B.49) 





Restricting the trace (B.l) to this background gives 


Tr e 


-sD 2 


M 


~ (4vr)- (d+1)/2 s [ dtd d xN [b 0 + Cl s 1/z a i( j + ...] 


(B.50) 


where the dots again indicate terms with more than two derivatives acting on the 
component fields. 

The computation of the heat-kernel coefficient C\ proceeds along the lines of App. 
B.2. The explicit form of the commutators entering into (B.25) is given by 


[ A t i A x ] — — (6a*a* + 2 ( did *) + 4a*(9*) A t , 

[At, [A*, A*]] =-8 a,a*A 2 , (B.51) 

[ Aj,, [ A t , A x ]] = - 16 aV didj A t + ... , 


where the dots again indicate terms which do not contribute to C\. Since multi¬ 
commutators of higher order automatically come with additional factors of a* or 
derivatives of a* the commutators (B.51) are the only terms that contribute to the 
heat-kernel coefficient c\. 

Abbreviating I = J dtd d xN a^a 1 the explicit evaluation of the traces (B.25) based 
on the commutators (B.51) yields 


11 IT (—) 

tl ~_—( rf+1 )/ 2 a-( d ~ 2 + z )/( 2z ) _ ^ 2z > T 

1_ ^ } z r(^A) ’ 


12 


T 2 ~ (47r) _(d+1)/2 s -(d-2+*)/(2*) / 2 r ( d 2z + X ) + ^ r ( d 2 z + x ) 


r (l) 


T(f) 


7 1 ~ ^ (AT T \-( d + 1 )/ 2 o-(. d ~ 2 + z )K 2z ) 7 ^ ( 2z ^ 2 ) r 

13 - 4 'A'v * ^ p ^d+2) 1 


Ti ~ - - (4vr)- (d+1)/2 s -(<i-2+d/(2d , / j ; 

3 


», . 7(^+4 
rm 


r (AA _l 11 

T b ~ (47r)"( d+1 ^ 2 s _ ( d_2+z )/( 22: ) 1 v 22 Z / 


r (I) 


Comparing the sum of these traces, 

13 


Tr e 


-sD 2 


f 4 vr'r (d+1)/2 s -( d ~ 2 + z )/( 2z ) - _ - r ( d 2z + il J 

6 d r(f) 


(B.52) 


(B.53) 


to the ansatz (B.50) finally gives the heat-kernel coefficient 


13 s - 1 r (^ +1) 
a d r (f) 


(B. 54 ) 


This result completes the derivation of the heat-kernel coefficients reported in Sect. 
3. 
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